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^ '. Abstract 



Inclusive cross sections for pion production in proton-proton col- 
. lisions are calculated for the first time based on unintegrated parton 

(gluon, quark, antiquark) distributions (uPDF). We use Kwiecinski 
uPDF's and phenomenological fragmentation functions from the lit- 
erature. In addition to the gg —* g diagram used recently for appli- 
cations at RHIC we include also gq —> q and qg — > q diagrams. We 
find that the new contributions are comparable to the purely gluonic 
one at midrapidities and dominate in the fragmentation region. The 
new mechanisms are responsible for ir + — tt~ asymmetry. We discuss 
^ ' how the asymmetry depends on ii? and pt- Inclusive distributions in 

xf (or rapidity) and transverse momentum for partons and pions are 
shown for illustration. In contrast to standard collinear approach in 
our approach the range of applicability can be extended towards much 
lower transverse momenta. 

PACS: 12.38.Bx, 13.85.Hd, 13.85.Ni 



1 Introduction 

The distributions of mesons at large transverse momenta in pp or pp collisions 
are usually calculated in the framework of perturbative QCD using collinear 
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factorization (see e.g. [H 121 EH E] ) • While the shape at transverse momenta 
larger than 2-4 GeV can be relatively well explained, the difference between 
the data and the lowest-order computation is quantified in terms of a so- 
called K-factor, independent quantity for each energy jlj. The K-factor is 
found to systematically decrease with growing energy [5, . In order to extend 
the calculation towards lower values of meson transverse momenta it was 
suggested to add an extra Gaussian distribution in transverse momentum 
jHlIZlE] 1 - In this approach the standard collinear integrals are replaced as: 



where the extra distribution function in transverse momentum is normalized 
to unity 



It is customary to use Gaussian distributions for g(k t ). It becomes clear that 
this procedure is effective in the following sense. The transverse momentum 
originates either from the nonperturbative "really internal" momentum dis- 
tributions of partons in nucleons (of the order of a fraction of GeV) and/or 
is generate dynamically as the inital state radiation process (of the order 
of GeV). In principle, the second component may depend on the values of 
longitudinal momentum fractions X\ and/or x%. The formalism used by us 
in the present paper will include both these effects separately and explicitly. 

The recent results from RHIC (see e.g. 13J have attracted a renewed 
interest in better understanding the dynamics of particle production, not only 
in nuclear collisions. Quite different approaches have been used to describe 
the particle spectra from the nuclear collisions [Hj. The model in Ref.|l()j 
with an educated guess for UGD describes surprisingly well the whole charged 
particle rapidity distribution by means of gluonic mechanisms only. Such 
a gluonic mechanism would lead to identical production of positively and 
negatively charged hadrons. The recent results of the BRAHMS experiment 
concerning heavy ion collisions [15J show that the n~ /n + and K~ / K + ratios 
differ from unity. This put into question the successful description of Ref . JO] • 
In the light of this experiment, it becomes obvious that the large rapidity 
regions have more complicated flavour structure. At lower energies these 
ratios are known to differ from unity drastically |28j . 

In Ref. [TTj one of us has calculated inclusive pion spectra in proton-proton 
collisions based on different models of unintegrated gluon distributions taken 
from the literature. In the present paper in addition to the gg — > g mecha- 

1 A similar procedure was used e.g. for prompt photon production [^j 






(2) 
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nism we include also q$g — > qj and gqf — > qj mechanisms and similar ones 
for antiquarks in order to obtain a fully consistent description. 

Many unintegrated gluon distributions in the literature are ad hoc parametriza- 
tions of different sets of experimental data rather than derived from QCD. 
An example of a more systematic approach, making use of familiar collinear 
distributions can be found in Ref.jTTJ. Recently Kwiecihski and collabora- 
tors [TBI [HI 1201 have shown how to solve the so-called CCFM equations 
by introducing unintegrated parton distributions in the space conjugated to 
the transverse momenta [IB]. We present first results for pion production 
based on unintegrated parton (gluon, quark, antiquark) distributions ob- 
tained by solving a set of coupled equations developed by Kwiecihski and 
collaborators. Recently these parton distributions were tested for inclusive 
gauge boson production in proton-antiproton collisions [21] and for charm- 
anticharm correlations in photoproduction [22]- While in the first process 
one tests mainly quark and antiquark distributions at scales /i 2 ~ M^,M|, 
in the second reaction one tests mainly gluon distributions at scales /x 2 ~ m 2 . 
In comparison to those reactions in the present application one tests both 
gluon as well as quark and antiquark distributions in a "more nonperturba- 
tive" region of smaller scales of the order down to fi ~ pt (parton) ~ 0.5 - 
1.0 GeV, which corresponds to pion transverse momenta pi(pion) ~ 0.25 - 
0.5 GeV. This is a region where perturbative and nonperturbative effects are 
believed to mix up and the application of the pQCD is doubtful. On the 
other hand, this is an interesting region of phase space responsible for the 
bulk of hadronic production. We shall discuss how far down to small pion 
transverse momenta one can apply the present approach. 

Some preliminary results of the present study were presented at a confer- 
ence [12]. 

2 Kwieciriski unintegrated parton distributions 

Kwiecihski has shown that the evolution equations for unintegrated parton 
distributions takes a particularly simple form in the variable conjugated to 
the parton transverse momentum. The two possible representations are in- 
terrelated via Fourier-Bessel transform 




(3) 
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/ fe (x,6,/i 2 ) 



dn t K t J (K t b)f k (x,K 2 t ^ 2 ) . 
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The index k above numerates either gluons (k=0), quarks (k> 0) or anti- 
quarks (k< 0). In the impact-parameter space the Kwiecihski equation takes 
the following simple form 

df N s{x,b,fi 2 ) a s {fi 2 ) H 



df s (x,b,fi 2 
d/j, 2 



27T/2 2 

- fNs(x,b,/j, 2 
a s (n 2 



dzP qq (z) 



0C f 

Q(z - x) J ((l - z)nb) f NS ( -, b, 



27T/2 2 



e{z-x)J {(l-z)fjb) 



P qq (z)f s (-,b,fj, 



+ P qg (z)f G [-,b,^ 



X 



[zP qq (z) + zP g q(z)} f S (x, b, fi 2 



df G (x,b,fi 2 ) a s (n 2 ) 



27T/X 2 



dz<Q(z-x) J ((l-z)fjb) 



Pgg(z)fs [-,b,H 



X 



+ P gg (z)f G [-,b,^ 



X 



[zPgg(z) + ZPgg(z)] f G (x, 6, // 



We have introduced here the short-hand notation 
/zvs — fu — fu, fd — fdi 

fs = fu + fu + fd + fd + fs + fs- 



(4) 



(5) 



The unintegrated parton distributions in the impact factor representation 
are related to the familiar collinear distributions as follows 



f k (x,b = 0,/i 2 



X . j 



(6) 



On the other hand, the transverse momentum uPDF's are related to the 
integrated distributions as 

"OO 

,2 r / ,2 



xp k {x,n 



dn t f k (x,K t ,fi ) . 



(7) 







While physically f k ( ) should be positive, there is no obvious reason 

for such a limitation for f k (x, b, fi 2 ). 

In the following we use leading-order parton distributions from ref.[29j 
as the initial condition for QCD evolution. The set of integro-differential 
equations in b-space was solved by the method based on the discretisation 
made with the help of the Chebyshev polynomials (see e.g. [20]) • Then the 
unintegrated parton distributions were put on a grid in x, b and fi 2 and the 
grid was used in practical applications for Chebyshev interpolation (see next 
section). 
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3 Inclusive cross sections for partons 



The approach proposed by Kwiecihski is very convenient to introduce the 
nonperturbative effects like internal (nonperturbative) transverse momentum 
distributions of partons in nucleons. It seems reasonable, at least in the first 
approximation, to include the nonperturbative effects in the factorizable way 



The form factor responsible for the nonperturbative effects must be normal- 
ized such that 



in order not to spoil the relation (JBJ). In the following, for simplicity, we use 
a flavour and ^-independent form factor 



which describes the nonperturbative effects. The Gaussian form factor in b 
means also a Gaussian initial momentum distribution oc exp(— k 2 6q) (Fourier 
transform of a Gaussian function is a Gaussian function). Gaussian form 
factor is often used to correct collinear pQCD calculations for the so-called 
internal momenta. Other functional forms in b are also possible. 

The gg —>■ g mechanism considered in the literature is not the only one 
possible. In FigQ we show two other leading order diagrams. They are 
potentially important in the so-called fragmentation region. The formulae 
for inclusive quark/antiquark distributions are similar to the formula for 
99 —> 9 JEj- The formulae for all the processes mentioned are listed below: 

for diagram A (gg— >g): 



f q (x,b, f , 2 ) = f^ CFM (x,b, f , 2 )-F^(b). 
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F NF (b = 0) = 1 



(9) 




(10) 



da A 16N C 1 



dyd 2 p t iV c 2 — lp 2 




(11) 



for diagram 



Bi (q/ g -»• q/): 

da Bl 16N C /4\ 1 



dyd 2 p t N 2 -l\9J p 2 t 



/ a s (n 2 ) / 9// i(xi,Ki,/i 2 ) / 9/2 (x 2 ,/t2,/i 2 ) 
/ J 

5 (2) (ki + k 2 - Pt) d 2 Kid 2 n 2 , 



(12) 
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for diagram B 2 (g q/ — > q/): 

da g2 16iV c /4\ 1_ 
dyrf 2 ^ ~ iV 2 - 1 \ 9j pf 

X] / «s( fi2 ) fg/l(Xl,l4,fJ?) f qf /2(x 2 ,K 2 2 ,fl 2 ) 

f J 

5 {2) (K 1 + K 2 -p t )d 2 K 1 d 2 K 2 - (13) 

These seemingly 4- dimensional integrals can be written as 2-dimensional in- 
tegrals after a siutable change of variables [TT] 

5 (2) (^l + £ 2 ~p t ) d 2 K Y d 2 K 2 = ! ... ^ . (14) 



The integrands of these "reduced" 2-dimensional integrals in q t = k[ — k 2 are 
generally smooth functions of qt and corresponding azimuthal angle <p qt . In 
the following we use two different prescriptions for the factorization scale /z 2 : 



• fi 2 = p 2 with freezing for p\ < fi 2 , 

• /i 2 = p 2 + 14. 

In Eqs. (jll|) . (|T2"j) and f)13|) the longitudinal momentum fractions 



«i/2 = 7= exp(±y) , (15) 



where m x is the effective mass of the parton. 

The sums in ()12j) and (jT5j) run over both quarks and antiquarks. The 
argument of the running coupling constant Q 2 above was not specified ex- 
plicitly yet. In principle, it can be p\ or a combination of pf, k 2 and k\- In 
the standard transverse momentum representation it is reasonable to assume 
Q 2 = min(p 2 , k 2 , k 2 ) (see e.g. |H]). In the region of very small p t usually 
p 2 < k\, k 2 and Q 2 = p 2 is a, good approximation. 

Assuming for simplicity that Q 2 = Q 2 (p 2 ) or p 2 (function of transverse 
momentum squared of the "produced" parton, or simply transverse momen- 
tum squared) and taking the following representation of the 8 function 

5 (2) (/«i +K 2 ~ Pt) = Jc^y J d2b ex P («i + ^2 - Ptjb , (16) 

the formulae (|12|) and ()13|) can be written in the equivalent way in terms 
of parton distributions in the space conjugated to the transverse momentum. 
The corresponding formulae read: 
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for diagram A: 

da A 16N r 1 



a«(Pt) 



dyd 2 p t N 2 — 1 pi 

f g /i(x 1 , b, /i 2 ) f g / 2 (x2, b, fx 2 )J (p t b) 2Mb , (17) 



for diagram Bi: 

da Bl 16iV c /4\ 1 

y /?//i( :r i> & >A* 2 ) f g /2(x 2 ,b,fi 2 )J (p t b) 2ixbdb , (18) 



/ 

for diagram B 2 : 



rfa 132 16iV^ /4\ 1 2 



dyd 2 pt iV c 2 - 1 V9/ pj 
^2 / /s/i^Ij^A* 2 ) fqf/2(x 2 ,b,fi 2 )J (p t b) 2-nbdb . (19) 

These are 1-dimensional integrals. The price one has to pay is that now the 
integrands are strongly oscillating functions of the impact factor, especially 
for large p t . The formulae (fTTj) . (fTHj) and (fT9"j) are very convenient to di- 
rectly use the solutions of the Kwiecihski equations discussed in the previous 
section. 

When extending running a s to the region of small scales we use a param- 
eter free analytic model from ref . j 



4 From partons to hadrons 

In Ref.|l()j it was assumed, based on the concept of local parton-hadron 
duality, that the rapidity distribution of particles is identical to the rapidity 
distribution of gluons. In the present approach we follow a different approach 
which makes use of phenomenological fragmentation functions (FF's). In the 
following we assume 9h = g . This is equivalent to r]h = r\ g = y g , where rjh 
and rjg are hadron and gluon pseudorapitity, respectively. Then 

y g = arsinh ( sinh y h J , (20) 
V Pt,h J 
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where the transverse mass mt,h — y m h + Pth- 111 or der to introduce phe- 
nomenological FF's one has to define a new kinematical variable. In accord 
with e + e~ and ep collisions we define a quantity z by the equation Eh = zE g . 
This leads to the relation 



Pt,h Jt \ 
Pt, g = J{m tih ,y h ) 

z 



(21) 



where the jacobian J(m t} h,yh) reads 
J(m t ,h,yh) = 



-1/2 



mi 



ml h cosh y h 



[22) 



Now we can write a given-type parton contribution to the single particle 
distribution in terms of a parton (gluon, quark, antiquark) distribution as 
follows 



dv p (r]h,Pt,h) 
drj h d 2 p t ,h 



dy p d 2 p t ,p 



Kvp ~ Vh) s 2 p t ,h 



zpu 
J 



dz D p ^ h (z,n 2 D ) 

d(T(y P ,Pt, P ) 
dy p d 2 p tt p 



(23) 



Please note that this is not an invariant cross section. The invariant cross 
section can be obtained via suitable variable transformation 



dy h d 2 p, 



9(vh,Pt,h) 



da p (y h ,p tA ) fd(y h ,p tjh )\ 1 dcr p (y h ,p tih ) 



dr]hd 2 p t ,i 



(24) 



where 



Vh 



ji *' + sinh Tj h + sinh r) h 

rt.h 



m l+Pjh 

p'i.h 



sinh 2 r/h — sinh r/h 



(25) 



Making use of the 6 function in ()23|) the inclusive distributions of hadrons 
(pions, kaons, etc.) are obtained through a convolution of inclusive distribu- 
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tions of partons and flavour- dependent fragmentation functions 

da(r] h ,p tA ) 



dr] h d 2 p t ,h 
dcr*(y g ,Pt, : 



+ D 1f^h{z,fi 2 D ) 
/=-3 



dy g d 2 p t ,g 



drr Bl 



dy qf d 2 p ti q 



+ D 1f-*h(z,H 2 D ) 



/=-3 



dy qf d 2 p t!q 



Pt,g = JPt,h/ z 



yq =r lh 
Vt,q = JPt,h/ z 



Pt,q = JPt,h/ z 



(26) 



One dimensional distributions of hadrons can be obtained through the inte- 
gration over the other variable. For example the pseudorapidity distribution 
is 



dojjih 
drj h 



d 2 Pt,i 



dcr(Vh,Pt,h) 



(27) 



dr]hd 2 Pt,h 

There are a few sets of fragmentation functions available in the literature 
(see e.g. El, E3)- 



5 Results 

As an illustration of the formalism, in the present paper, we shall show results 
for energies adequate for CERN SPS i.e. for the energies at which the missing 
mechanisms should play an important role. 

Before we adress the distributions of pions we wish to discuss the inclusive 
spectra of "produced" partons. 

5.1 Parton distributions 

In the familiar collinear approach the contributions of diagrams involving 
quarks and antiquarks are not negligible even at large energies. A nice, quan- 
titative discussion of this issue can be found in [3]. In this section we shall 
make a similar analysis in our approach based on CCFM uPDF's. In FigO 
we display the contributions of partons (gluons, quarks, and antiquarks) for 
all diagrams of Fig^for the center-of-mass energy W = 17.3 GeV, i.e. energy 
of recent experiments of the NA49 collaboration at CERN as a function of 
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jet (minijet) Xp. The parton rapidity is related to the parton xp as follows 



V , 7 , (28) 

yJml + x 2 F l+p 2 t -x F f J 

where the parton mass here is the same as in Eq.(|15|). The corresponding 
cross section is obtained by integration over parton transverse momenta in 
the interval 0.2 GeV < pt < 4 GeV. The gg — > g contribution, claimed 
to be the dominant contribution at RHIC ^0], is somewhat smaller than 
the contribution of diagrams B\ and B 2 . The contribution of diagram B\ 
(dashed line) dominates at negative Feynman-x^, while the contribution of 
diagram B 2 (dotted line) at positive Feynman-x^. By symmetry require- 
ments da B2 1 'dxp{xp) = da Bl I dxp{— xp). 

In order to understand the intriguing asymmetry in xp of contributions 
of diagrams B\ and B2, in FigOlwe present a further decomposition of the B\ 
contribution into sea-glue, valence-glue sub contributions and correspondingly 
of the B 2 contribution into glue-sea, glue-valence subcontributions. This de- 
composition shows that the sea-glue and glue-sea contributions are of similar 
size as the valence-glue and glue-valence, respectively. It is interesting to 
note a different a^-asymmetry of the contributions involving sea and valence 
quarks. This decomposition explains the shift of maxima of contributions 
corresponding to diagram B\ and B 2 seen in FigJ21 

For completeness in FigEJwe show transverse momentum distribution of 
"produced" gluons corresponding to diagram A and of quarks and antiquarks 
corresponding to diagrams B\ and B 2 . In this calculation we have integrated 
over parton xp (or parton rapidity). For the two contributions one obtains a 
rather similar functional behaviour. It is worth noticing that in contrast to 
the standard collinear case in our approach the partonic cross section is fully 
integrable. This seems pro missing in extending the region of applicability of 
the "perturbative" 2 QCD towards smaller transverse momenta of hadrons. 

The rise of the cross section above p t ~ 0.5 GeV is due to a rapid increase 
of gluonic radiation above p t = /io, where /i 2 , is the minimal factorization scale 
for the GRV parton distributions. On the other hand, the rise of the cross 
section towards pt < fJ>o is an artifact of freezing factorization scale below //q 
in the conjuction of the singular behaviour of the denominator in the parton 
cross section formulae. To elucidate this issue somewhat better in FigElwe 
present also the results with the second prescription for the factorization 

2 Our approach is not completely perturbative. Some nonperturbative effects are con- 
tained in the nonperturbative form factor, way of freezing a s or the scale of parton distri- 
butions. 
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scale with the minimal shift of the factorization scale (dotted line) and the 
result with an extra substitution of 4^ by p i+ m s (dashed line), where m 2 d is 
a new nonperturbative parameter of our model specified in the figure. The 
second low-scale prescription lead to a smooth behavior of the partonic cross 
section at low transverse momenta. The consequences and the sensitivity to 
the extra prescription on the pion production at low transverse momenta of 
hadrons will be discussed in the next section. 

5.2 Pion distributions 

The passage from parton distributions to hadron distributions is the next step 
of our analysis. In the present approach we shall use the approach based on 
phenomenological fragmentation functions known from other processes. In 
the literature it is mostly e + e~ — > hadrons reactions which are used for 
extraction of phenomenological fragmentation functions. 

There are a few sets of fragmentation functions in the literature. In the 
present calculation we shall use the fragmentation function of Kretzer [25] . 
There are two advantages of this particular set of fragmentation functions 
over other known from the literature. Firstly, the Kretzer fragmentation 
functions are available and reasonable even for very low scales p t ~ 1 GeV, or 
even less, i.e. in the region of our interest. Secondly, more attention in their 
construction was paid to the flavour structure than in any other approach 
in the literature. A good control of the flavour structure is essential when 
discussing and comparing contributions of diagram A, Bi and B 2 to the pion 
momentum distributions. 

In FiglHl we compare our model invariant cross sections for pp — > tt + 
(left panel) and pp — > n~ (right panel) as a function of pion transverse 
momentum at W = 27.4 GeV for different values of the parameter bo of our 
Gaussian nonperturbative form factor. In principle, our result should not 
exceed experimental data especially in the perturbative regime of p t > 2 
GeV where the perturbative 2^2 parton subprocesses are crucial. This 
limits the value of the nonperturbative form factor to &o > 0.5 GeV -1 . 

As can be seen from the figure the agreement of our model with exper- 
imental data is not perfect. There can be a few reasons for this fact. In 
order to understand the problem somewhat better let us concentrate on a 
small transverse momentum region p t ^ < 2 GeV. As an example in Fig|7| 
we present our results for b = 1 GeV -1 . We observe a deficit at p t > 0.5 
GeV and a strong excess at p t < 0.3 GeV. It is particularly intriguing which 
effect stands behind the huge excess at very small transverse momenta. In 
addition, in Fig|7|we present contributions of different disjoint and comple- 
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mentary regions of variable z, as indicated in the figure. The figure shows 
that the rapid increase below p t ~ 0.3 GeV is caused exclusively by small z < 
0.2 i.e. must be traced back to the phenomenological fragmentation functions 
used. The Kretzer fragmentation functions as calculated from the publicly 
available computer code [2H] are restricted to factorization scales larger than 
1 GeV 2 . Therefore in our calculation we were forced to freeze the fragmenta- 
tion functions below — Pt < 1 GeV. The fragmentation functions at this 
small scale are shown by the solid line in FiglHlfor illustration. A clear rise of 
the fragmentation functions at small z can be observed, which in conjunction 
with the previous observation must be responsible for the excess of pions at 
very small p t ,h- The lower limit for the scale of about 1 GeV 2 is a recomenda- 
tion based on analysis of the e + e~ data rather than a rigorous applicability 
limit. In fact, the initial evolution scale in [23] is /j,q = 0.26 GeV 2 . Let us try 
to use the Kretzer fragmentation functions at this even lower scale. In FigGO 
we show (dashed line) the gluon and quark fragmentation functions. While 
for g — > 7T fragmentation the small-z rise in "z D(z)" completely dissapears 
and a minimum at z = is obtained (valence-like character), for the quark 
fragmentation only a saturation of "z D(z)" can be observed. In FigE] we 
show the result for pion production obtained with the fragmentation function 
at the lowest factorization scale /Mq = 0.26 GeV 2 . Clearly the description of 
the data is now better. It is somewhat amusing that even the description 
above p t ,h = 0.5 GeV is now much better. It becomes obvious that the choice 
of the factorization scale is very important in this context. 

The description below p t ,h — 0.2 GeV is still rather poor. This is partially 
caused by the behaviour of the quark fragmentation functions at small z (see 
FiglHJ). Other effects are also possible. For example, we have used a simple 
one-parameter Gaussian form factor. In principle, such form factor may have 
more complicated functional dependence and can depend not only on the im- 
pact parameter, but also on energy and/or x\ and x 2 - Furthermore the whole 
formalism with z-dependent fragmentation functions must ultimately break 
at very low transverse momenta, in particular at small xf where resonance 
decays must be treated explicitly. Having this in view our final result seems 
promissing for further improvements in the future. 

Inclusion of diagrams B\ and B 2 in conjunction with the flavour depen- 
dent fragmentation functions lead to the n + — 7r~ asymmetry. In FigHT)! we 
show the asymmetry as the function of pion transverse momentum. The 
asymmetry is well described by our model, in contrast to individual distri- 
butions. This seems to suggest the right relative contributions of diagram 
A, Bi and B 2 . The asymmetry depends only weakly on the value of the 
parameter bo of the Gaussian nonperturbative form factor. 

At SPS energies it is often customary to use xf rather than rapidity. The 
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hadron Feynman variable xp,h is related to hadron pseudorapidity rjh as 



In Fig^^we compare the xp distributions obtained in our approach for 6 
= 1 GeV 2 with those from some popular programs available on the market. 
The thin and thick (dashed and dotted) lines correspond to the different 
treatment of the scale in the fragmentation functions as described above (see 
also the figure caption). Clearly some mechanisms at large xf are missing 
in our approach. One of such mechanisms is the so-called pion stripping 
discussed e.g. in Ref.[3U]. It was not our goal to include all of these effects 
and describe the experimental data. Our intention here was to identify the 
dominant mechanisms in the language of unintegrated parton distributions 
adequate at intermediate and low transverse momenta. 

For completeness in Fig^] we present a decomposition of the pion xp 
distribution into the components corresponding to diagrams shown in Figd 
The purely gluonic contribution (diagram A) is comparable to the other two 
contributions at xp ~ 0, but at \xp\ > 0.5 the contributions of diagram 
Bi and B 2 are significantly larger. Similarly as for parton distributions the 
diagrams B\ and B2 contribute at somewhat larger \xp\ than the diagram A. 
Some differences for 7r + and tt~ distributions are visible in the figure. The 
mechanisms underlying diagrams Bi and B 2 are important contributors to 
understand the 7r + — 7r~ asymmetry at small and intermediate values of xp. 
We think, however, that the dominant mechanism of ir + and n~ asymmetry 
in very forward (xf > 0.5) and very backward (xf < -0.5) hemispheres is 
due to the pion stripping mechanism (for a reference see e.g. |3T)]). 

6 Conclusions 

The approach based on unintegrated gluon distributions was applied recently 
to describe particle momentum distributions at RHIC for nucleus-nucleus JTU] 
and proton-proton collisions JT] . We propose new mechanisms, neglected so 
far in the literature, which involve also quark/antiquark degrees of freedom 




(29) 



and the jacobian of transformation expressed by xp,h and pt,h is 




(30) 
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and are based on (anti)quark-gluon and gluon-(anti) quark fusion processes 
followed by the subsequent fragmentation. These missing mechanisms have 
been estimated in the approach based on unintegrated parton (gluon, quark, 
antiquark) distributions originating from the solution of a set of coupled 
equations proposed recently by Kwiecihski and coworkers. The formalism 
proposed recently by Kwiecihski is very useful to obtain not only gluon unin- 
tegrated distributions but also their counterparts for quarks and antiquarks. 

In the present paper we have concentrated on low energies W ~ 20 GeV, 
relevant for SPS experiments. By freezing relevant scales for small parton 
transverse momenta we achieve a reasonable description of the data down to 
pion transverse momenta ~ 0.5 GeV. The missing terms lead to an asymme- 
try in the production of the 7r + and ir~ mesons. While at lower energies of the 
order of 20 - 50 GeV such an asymmetry has been observed experimentally, 
it was not yet studied at RHIC. The BRAHMS collaboration at RHIC has a 
potential to study the asymmetries. Such asymmetries may be also useful to 
pin down the non-gluonic mechanisms as those encoded in diagrams B\ and 
B 2 . 

The non-ideal agreement with the experimental data at low p t can be due 
to approximate treatment of nonperturbative effects embodied in the form 
factors, as well as in the treatment of hadronization effects with the help of 
scale-dependent one-parameter fragmentation functions. Both these effects 
require further detailed studies which go beyond the scope of the present 
paper. 
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Figure 1: Leading-order diagrams for inclusive parton production 
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Figure 2: Distribution of virtually "produced" partons as a function of xp 
for W = 17.3 GeV and bo = 1 GeV^ 1 of the Gaussian form factor. In this 
calculation 0.2 GeV < p t < 4 GeV. Contribution of diagram A is shown by 
thin solid line, the contribution of diagram B\ by dashed line and contribution 
of diagram B2 by dotted line, and the sum of all processes by thick solid line. 
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Figure 3: Decomposition of diagram B\ cross section into glue-sea and 
glue-valence components and of diagram B 2 cross section into sea-glue and 
valence-glue components as a function of xf- 
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Figure 4: Transverse momentum distribution of partons from diagram A (left 
panel), and B\ + B 2 (right panel) for W = 27.4 GeV and different values of 
the Gaussian form factor. The cross sections presented were obtained through 
integration over -1 < xp < 1. 
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Figure 5: Transverse momentum distribution of partons from diagram A (left 
panel), and B\ + B 2 (right panel) for W = 27.4 GeV and -1 < xf < 1- In 
this calculation bo = 1 GeV^ 1 . The solid line is the same as in the previous 
figure, i.e. freezing prescription for the factorization scale was used. For 
comparison the dotted line includes the shift prescription for the factoriza- 
tion scale. Finally the dashed line includes shift of factorization scale and 
modification of denominator as described in the text. 
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Figure 6: Invariant cross section as a function of transverse momentum of 
7T + (left panel) and tt~ (right panel) for rj = 0, W = 27.4 GeV and different 
values of parameter in the Gaussian form factor. Experimental data for W = 
23, 31 GeV from \27\j and for W = 27.4 GeV [28] are shown for comparison. 
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Figure 7: Invariant cross section as a function of transverse momentum of 
7r + (left panel) and n~ (right panel) for r) — 0, W = 27.4 GeV. In this 
calculation b = 1 GeV -1 . 
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Figure 8: The basic Kretzer fragmentation functions at low factorization 
scales. 
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Figure 9: Invariant cross section as a function of transverse momentum of 
7T + (left panel) and n~ (right panel) for rj — 0, W = 27.4 GeV obtained with 
the Kretzer fragmentation functions at /!q = 0.26 GeV 2 . In this calculation 
b = 1 GeV- 1 . 
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Figure 10: Ratio of cross sections from Fig\$ as a function of pion transverse 
momentum for W = 27.4 GeV. 
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Figure 11: Our result against the UrQMD (red triangles) and HIJING (blue 
squares) models for SPS energy W = 17.3 GeV, for positive (dashed) and 
negative (dotted) pions. Thin lines correspond to FigFft and thick lines to 
Fig\B 
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Figure 12: Contribution of diagrams A (thin solid), B\ (dashed) and B2 (dot- 
ted) to charged-pion xf distributions for W = 17.3 GeV. In this calculation 
b = 1.0 GeV- 1 and 0.2 GeV < p t;h < 2 GeV. 
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